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A  HYBRID  MODE  AND  A  CLASSIFICATION 
OF  BEAM  PLASMA  INSTABILITIES 

I.  Introduction 

There  is  considerable  renewed  interest  in  the  old  subject  of  the  electron 
stream  interactions  with  plasmas.  While  the  classical  two  stream  instability 
is  quite  well-understood,  much  of  the  current  effort  has  been  placed  on 
intense  relativistic  electron  beams.  Examples  of  areas  of  interest  where  such 
a  study  is  warranted  include;  electron  beam  heating  of  a  dense  plasma  for 
inertial  confinement  fusion  driver;^  laser-plasma  beat  wave  acceleration  of 
particles  to  ultra  high  energies; electron  beam  propagation  under  sub- 

ll  c  c. 

atmospheric  conditions"^  and  in  the  case  of  a  modified  betatron  accelerator, 
the  proposed  method^  of  injection  and  extraction  of  the  electron  beam  by 
creating  a  plasma  channel;  etc. 

The  diverse  areas  mentioned  above  encompass  a  multi-dimensional  parameter 

space  which  spans  over  many  orders  of  magnitude.  Depending  on  the  objective 

of  an  investigation,  a  theoretical  analysis  of  beam-plasma  interaction  is' 

necessarily  limited  in  scope,  as  the  electromagnetic  effects,  the  relativistic 

effects,  and  the  coupling  between  the  perpendicular  and  parallel  motions, 

among  different  species  of  particles,  compete  with  each  other.  In  this  paper, 

ft"  1  2 

we  attempt  to  classify  various  beam-plasma  instabilities  according  to  the 
electron  beam  energy  and  the  electron  beam  density  (relative  to  the  background 
plasma  density).  Domains  of  instability  mechanisms  are  determined,  within 
each  of  which  the  instability  growth  rate  and  the  minimum  axial  wavelengths 
for  the  instabilities  are  displayed.  Such  a  classification  then  provides  an 
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immediate  determination  of  the  nature  and  the  importance  of  a  particular 
instability,  once  the  beam  energy,  the  beam  density,  and  the  background  plasma 
density  are  given.  This  would  allow  more  ready  refinement  of  the  theory,  if 
necessary,  as  the  dominant  physical  mechanisms  are  isolated  in  such  a 
classification. 

During  our  survey  of  the  parameter  space,  we  identified  an  instability, 
which  we  termed  a  hybrid  instability.  This  instability  is  important  for  the 
propagation  of  a  long,  thin  relativistic  electron  beam  into  a  tenuous 
background  of  ions.  The  phase  (and  group)  velocity  of  the  wave-  along  the 
beam,  in  the  laboratory  frame,  is  1/2  of  the  speed  of  light.  The  growth  rate 
is  proportional  to  the  quarter  powers  of  the  electron  density  and  of  the  ion 
density,  and  inversely  to  the  quarter  powers  of  beam  energy  and  of  the  ion 
mass . 

The  hybrid  instability  could  be  important  in  the  "ion  focus  regime", 
where  an  intense  electron  beam  propagates  in  a  plasma  channel,  which  is  either 
preformed  or  self -generated.  If  the  plasma  channel  is  tenuous,  the  beam 
electrons  electrostatically  eject  the  plasma  electrons  from  the  beam  path, 
leaving  behind  the  plasma  ions  which  then  guide  the  beam  electrostatically. 

An  example  where  this  hybrid  instability  is  possibly  a  dominant  one,  is  when 
a  preformed  plasma  channel  is  generated  to  extract  an  electron  beam  in  a 
modified  betatron  configuration,  after  the  electron  beam  is  fully  accelerated. 
(See  examples  below). 

In  this  paper,  we  shall  first  present  a  simple  physical  argument  which 
leads  to  the  dispersion  relationship  of  the  hybrid  mode  for  a  slab  model.  It 
turns  out  that  the  simple  dispersion  relationship  thus  derived  is  an  excellent 
approximation  to  the  full  dispersion  relationship,^  ^  which  incorporates  all 
electromagnetic  and  relativistic  effects.  We  find  that  the  dispersion 


relationship  is  also  applicable  to  a  cylindrical  beam  of  small,  but  finite 
cross  sections  [see  Appendix]. 

To  place  the  hybrid  mode  in  a  proper  perspective,  in  Section  III  we 
provide  the  classification  scheme,  based  on  the  beam  energy  and  beam  density. 
This  classification  identifies  the  domains  for  the  hybrid  mode,  the  Weibel 
mode,^^  and  the  classical  two  stream  instabilities.  In  that  section,  we  also 
furnish  a  simple  derivation  of  the  Weibel  mode  for  a  relativistic  electron 
beam.  The  growth  rates  and  the  minimum  axial  wavelengths  for  the  various 
instabilities  are  catalogued  for  easy  reference.  Several  examples  of  current 
interest  are  given  in  Section  IV  to  demonstrate  the  ready  utilization  of  the 
classification  scheme  established  in  Section  III. 


The  basic  assumptions  are  now  stated.  We  assume  that  the  external 
magnetic  field  is  unimportant.  This  is  certainly  true  in  the  ion  focus  regime 
mentioned  above.  The  plasma  and  the  beam  are  assumed  collisionless.  The 
plasma  is  at  rest,  and  both  the  plasma  and  the  beam  are  cold.  It  is  assumed 
that  there  is  no  plasma  return  current.  One  may  argue  that  such  an  assumption 
is  valid  in  the  ion  focus  regime  where  the  background  plasma  density  is  low. 
The  classification  is  based  on  a  model  of  infinite  uniform  beam  plasma 


system.  The  latter  assumption  is  relaxed  in  the  Appendix.  In  fact,  if 
2  2  2 

k  -  1 /r.  >>  k  where  k  is  the  perpendicular  wave  number,  k  is  the  axial  wave 
1  ^  1 

number  and  r^^  is  the  electron  beam  radius,  the  slab  model  provides  an  adequate 


description  for  global  modes  of  a  sharp  boundary  model  if  is  taken  as  a 

suitably  quantized  number  greater  than  unity  in  order  of  magnitude. ^ ^  It 


is  within  this  context  that  our  classification  of  instabilities  is  made.  We 


shall  further  discuss  these  basic  assumptions,  and  the  limitations  they 


impose,  in  the  last  section 


II .  The  Hybrid  Mode 


The  dispersion  relationship  for  the  hybrid  mode,  together  with  the 
physical  picture,  is  now  presented.  For  simplicity,  we  consider  a  one¬ 
dimensional  model.  The  approach  is  very  similar  to  the  one  used  to  describe 
simple  plasma  oscillations. 

Consider  an  infinitely  uniform  electron  beam  propagating  in  the  z- 
direction  with  constant  velocity  v^.  The  beam  is  highly  relativistic,  i.e., 

y  »  1  (1  ) 

2  2  1/2 

where  73(1-  v^/c  )  is  the  relativistic  mass  factor  and  c  is  the  speed  of 

light.  The  beam  is  cold  and  is  embedded  in  a  uniform  background  of  ions  of  , 

( 

considerably  lower  density.  The  ions,  of  mass  m^ ,  are  also  cold  and  are 
motionless  in  equilibrium.  We  assume  that  there  is  a  force  balance  for  the 
ions  and  for  the  electrons  in  equilibrium.  For  the  electrons  the  almost 
perfect  cancellation  between  the  self  electric  field  and  the  self  magnetic 
field,  (for  Y  >>  1 )  together  with  the  weak  ionic  background,  would  make  a 
perfect  equilibrium  of  cold  electronic  flow  plausible.  For  the  ions,  because 
of  their  mass,  we  simply  ignore  whatever  ionic  motions  which  may  have  resulted 
from  the  equilibrium  electric  field.  If  the  characteristic  time  of  the 
instability  is  short  compared  with  the  free-fall  time  of  the  ions,  one  may 
then  ignore  the  external  forces  on  the  ions  in  the  equilibrium  states,  to  the 
lowest  approximation.  This  in  fact  can  be  shown  to  be  the  case  through  a 
Vlasov  description  of  the  ionic  responses.  The  effects  of  these  equilibrium 
ion  motions  may  then  be  evaluated  in  a  qualitative  manner  by,  for  example, 
pretending  that  there  is  a  "temperature”  of  the  ions  for  the  equilibrium. 

Here,  we  just  assume,  for  convenience,  that  such  temperature  is  zero.  Simply 
stated,  we  assume  that  the  unperturbed  equilibrium  fields  are  unimportant. 
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This  is  certainly  true  in  the  hypothetical  situation  of  a  completely 


neutralized  beam  and  plasma.  The  above,  then,  constitutes  the  working 
hypothesis. 

To  consider  the  small  signal  stability  of  such  a  beam-plasma  system,  we 
focus  mainly  on  the  transverse  motion  of  the  beam  since  a  beam  electron  is 
more  responsive  to  a  transverse  force  than  to  a  longitudinal  force  when  the 
beam  is  highly  relativistic.  Let  be  the  electron  displacement,  from  its 
equilibrium  position,  in  the  transverse  direction,  and  be  the  ion 
displacement,  also  in  the  transverse  direction  [Fig.  1].  These  displacements 
then  generate  AC  electric  fields  E^g  and  and  AC  magnetic  fields  B^^g  and 

.  Here  the  subscripts  e  and  i  are  used  to  denote  that  part  of  electric 
field  (or  magnetic  field)  which  is  due  to  5^  and  respectively.  Thus 


Eie  -  4,r,e|nj^C^ 


=li  = 


(2) 


B.  -  -E,  V  /c, 
le  le  o 


^i 


(3) 


where  n^j  and  np  are,  respectively,  the  beam  electron  and  plasma  ion  density  in 
equilibrium  and  e  is  the  electron  charge.  Here  we  assume  that  the  ions  are 
singly  ionized. 

Equation  (2)  is  obvious.  B^^g  as  given  by  Eq.  (3)  may  easily  be  obtained 
by  visualizing  the  dipole  current  layers  labeled  0  and  E  in  Fig.  1  .  Since 
there  is  no  unperturbed  ion  motion,  B^^^^  =  0,  as  shown  in  Eq.  (3). 

The  linearized  equations  of  motions  for  the  electrons  and  ions  then  read 


le 


li 


(4) 
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"  leiCE^^  -  E^.]. 


Since  we  assume  that  there  is  no  ionic  motion  in  equilibrium,  the  Lorentz 

force  is  absent  for  the  ions  [cf.  Eq.  (5)].  Further  simplification  is 

possible  by  using  Eq.  (1)  and  by  the  assumption  that  the  background  plasma  is 

tenuous.  By  Eq.  (1),  E,  +  (v  /c)  x  B,  =0  since  the  self  electric  field  and 

le  o  le 

the  self  magnetic  field  generated  by  the  electronic  dipole  layer  cancelled  to 
order  1/Y^,  and  Y  is  large  [cf.  Fig.  1;  also,  Eqs.  (2),  (3)].  The  latter 
assumption  permits  us  to  ignore  E^^j^  in  comparison  with  E^g  in  Eq.  (5).  Thus, 
Eqs.  (i<)  ana  (5)  are  simplified  to  read 


Y  m  £  *  -  leiE,.  =  ‘^irtei  n  5- 
o  e^e  '  I  li  '  '  p^i 


i.^i-  -  lejE^g  =  *  i)ir,e|  n.^5, 


For  a  perturbation  with  dependence  exp  (iut  -  ikz),  =  -  (co  -  kv^)  and 
2 

=  -  0)  Using  these  expressions  in  Eqs.  (6)  and  (7)  and  multiplying  the 

resultant  equations,  we  obtain  the  dispersion  relationship  for  the  hybrid  mode 


2  2  2  2 
0)  (o)  -  kv  )  =  0)  (jo, 

o  P  t)J_ 


where 


2  2 
u)  a  4ire  n  /m. 
P  P  1 


»  ilire^n. /m  Y. 
0j_  b  e 
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The  instability  growth  rate  may  readily  be  obtained  from  a  square  root  of  Eq. 
(8): 


(u  =  TT  kv  ±  i  (o) 
2  0 


/  , 2  2,1/2 


(11) 


This  dispersion  relationship  has  been  shown  to  be  an  excellent  approximation 

to  that  for  a  thin  beam  which  is  highly  relativistic  (B  1).  The  latter  is 

given  by  Eq.  (All)  of  the  Appendix.  The  solution  to  Eq.  (All)  is  slightly 

more  general  than  (11)  and  is  sketched  in  Fig.  2,  where  the  real  part  and 

2  2  2  2  2 

the  imaginary  part  ui.  are  shown.  In  Fig.  2,  n  =  =  o).  /"Y  o)  .  Equation 

1  OF  [j21  p  bj_  p 

(All)  reduces  to  Eq.  (11)  in  the  limit  S  =  1. 

Several  points  for  the  hybrid  mode  are  noteworthy:  (a)  The  phase  and 
group  velocity  of  the  unstable  waves  equal  to  v^/2  =■  c/2,  (b)  The  instability 

is  absent  if  the  axial  wavelength  A  satisfies  the  inequality 


2tr  ^  ire 

/o)  0), 

P  Dl 


(12) 


Here,  we  have  used  a  subscript  H  to  denote  the  critical  axial  wavelength  for 
the  hybrid  mode,  (c)  The  instability  growth  rate,  for  k  <<  k^^,  is  given  by 


iH 


-  /oj  OJ, 

P  bl 


(13) 


which  is  independent  of  k^  and  of  k.  (d)  We  may  more  precisely  state  what  we 
meant  by  a  "tenuous"  ionic  background.  For  <<  as  we  have  assumed 

in  deriving  (7)  from  (5),  we  must  have  >>  [cf.  (2)].  This 

inequality  becomes 
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2  2 

<*>K  »  1 

bj_  P 


(1'4) 


If  we  use  in  (7).  (e)  While  the  dispersion  relationship  (11)  is 

independent  of  k^,  the  transverse  displacement  is  extremely  important.  An 

analogous  example  would  just  be  the  simple  plasma  dispersion  relation 

(u^  -  u)^.  which  is  independent  of  k  (but  displacements  of  particles  parallel 
P 

to  k  are  important). 

At  the  risk  of  some  redundancy,  we  paraohrase  the  above  derivation  for 

the  hybrid  mode  by  a  physical  description  of  the  instability  mechaniom,  as 

depicted  in  the  sequence  of  events  in  Fig.  3*  Suppose  that  a  sheet  of 

electrons,  whose  unperturbed  position  is  on  the  sheet  labeled  0  as  shown  in 

Fig.  3(a),  is  displaced  to  position  labeled  E,  at  a  distance  Cg  from  0. 

Assume  also  that  a  sheet  of  ions  is  displaced  to  level  I,  by  a  distance 

from  0.  Let  C-  >  C  =^3  shown  in  Fig.  3(a)-  If  Lhe  electron  beam  is  highly 
i  16 

relativistic,  the  force  on  the  electron  sheet  E  is  due  only  to  the  electric 
field  Ej^j^  generated  by  the  ion  sheet  I,  since  the  self  electric  field  Ej^g  and 
the  self  magnetic  field  B^^g,  which  are  generated  by  the  dipole  sheets  E  and  0, 
are  canceled  to  order  1/Y^.  Thus  sheet  E  is  accelerated  toward  sheet  I,  due 
only  to  the  ionic  electric  field. 

Now  consider  the  forces  on  the  ion  sheet  I  in  Fig.  3(a).  Sheet  I  is  not 
acted  on  by  the  electrons  at  all,  since  the  electric  field  E^^g  and  magnetic 
field  Bj^g  (generated  by  E^)  are  confined  within  the  region  between  E  and  0. 
(One  may  visualize  this  by  imagining  sheet  E  and  0  as  capacitor  plates,  which 
generate  Ej^g,  and  also  as  a  dipole  sheet  current,  which  generated  B^^g).  Thus 
sheet  I  is  subject  only  to  its  own  restoring  force  However,  since  the 

ion  background  is  tenuous,  E^^j^  is  small.  In  other  words,  the  ion  sheet  I  in 
Fig.  3(a)  is  not  subject  to  any  perturbed  force,  to  the  lowest  order.  Since  E 
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is  accelerated  toward  I,  as  explained  in  the  previous  paragraph,  E  overtakes  I 

at  a  later  time,  as  shown  in  Fig.  3(b). 

In  Fig.  3(b),  >  5^.  Now  sheet  E  is  not  subject  to  any  perturbed 

force,  to  the  order  we  are  considering,  since  E^^g  and  canceled  to  order 
2 

1/Y  ,  and  E^j^  is  confined  only  within  I  and  0.  Thus  once  E  surpasses  I,  as  in 

Fig.  3(b),  E  moves  with  a  constant  velocity.  However,  I  is  accelerated  toward 

E,  because  of  E^^g  (again,  Ej^j^  is  small  because  of  the  low  background  ion 

density).  Thus,  at  a  later  time  I  overtakes  E,  as  shown  in  Fig.  3(c),  and  the 

whole  sequence  repeats.  Thus  the  electron  sheet  E  and  the  ion  sheet  I 

gradually  deviate  from  their  original  position  0,  and  this  is  essentially  the 

basic  mechanism  which  drives  the  hybrid  mode. 

The  above  physical  description  is  in  some  respect  similar  to  the  one 

given  for  simple  plasma  oscillations  (u  =  ujp).  For  example,  the  argument  does 

not  require  a  detailed  specification  of  the  wave  numbers,  of  the  beam  and 

plasma  dimensions  and  of  the  profiles.  One  may  then  be  tempted  to  conclude 

that  the  hybrid  instability  is  a  fundamental  one  as  long  as  >>  1  and 
2  2 

>>  1 .  In  fact,  we  shall  show  in  the  Appendix  that  essentiaxly  the  same 
dispersion  relationship  (8),  or  (11),  is  obtained  for  a  cylindrical  beam  with 
a  sharp  boundary  profile. 

We  postpone  to  the  next  two  sections  for  a  discussion  of  the  relationship 
of  this  hybrid  mode  and  other  beam-plasma  instabilities. 

III.  A  Classification  of  Beam  Plasma  Instabilities 

To  place  the  hybrid  mode  in  the  "general  picture"  of  beam  plasma 
interactions,  we  have  conducted  a  limited  search  of  the  literature  in  which 
the  relativistic  and  electromagnetic  effects  have  been  treated.  We  found  that 
the  recent  work  of  Tajiraa,^^  Jones, and  Shoucri^^  are  the  closest  to  the 


present  one.  The  studies  of  Jones  and  Tajima  are  more  general,  as  they 
included  bounded  systems  in  cylindrical  geometries.  Their  results  firmly 
establish  the  validity  of  the  use  of  a  slab  model,  if  the  beam  radius  r^j  (or 
thickness  t)  is  sufficiently  small,  ^»nd  if  the  perpendicular  wave  number  in 
the  slab  model  is  to  be  replaced  by  where  p  is  a  quantized  number  of 

order  unity,  (or  larger  for  higher  radial  modes).  Shoucri  considered  only  a 
slab  geomtry,  but  he  brought  out  an  interesting  point  regarding  the  coupling 
of  electromagnetic  and  relativistic  effects.  We  postpone  to  Section  V  for  a 
brief  discussion  of  other  literature. 

The  three  authors  mentioned  in  the  previous  paragraph  produced,  among 
other  things,  the  same  dispersion  relationship.  For  example,  the  dispersion 
relationship  given  in  Eq.  (20)  of  Shoucri  is  equivalent  to  the  dispersion 
relationship  (40)  of  Tajima,^*^  which  is  in  turn  equivalent  to  that  given  in  p. 
1931  of  Jones, who  also  noted  such  a  correspondence.  This  dispersion 
relationship,  for  the  slab  model,  reads 

B  E  =  C  (15) 

where 

B  =  1-  u)p/uj^  -  n6) 

t7  1  .22,2  ,2  2,  2  ,2^  2.,  2 

E-1-kc/u)  -  kc/uj  “(a>  '♦■oj.  )/a)  (17) 

I  P 

^  .2  2  2  ,  2  2.  ,  2_2 

C  -  k  c  u  (m.  n  ~  «.  )/u)  0  •  (18) 

j_  p  bll 

Here,  B  represents  the  classical  beam  plasma  interaction,  E  the 
electromagnetic  mode,  and  C  the  coupling  of  the  two.  In  Eqs.  (16)-(18),  u 


is  the  frequency,  k  is  the  axial  wave  number,  is  the  perpendicular  wave 

number,  ft  -  u  -  kv^,  is  the  plasma  frequency  of  the  background, 

2  2 

=•  4ire  nj^/m^Y,  and  C  -  0  if 

k  »  0,  or  if  0)  =  0,  or  in  the  non-relativistic  limit,  in  which  case 

1  P 

Since  the  dispersion  relation  (15)  is  derived  under  the 

assumption  that  the  unperturbed  force  on  each  of  the  particle  species  is 

negligible,  hereafter  oi  may  stand  for  ui  .  or  oj  depending  on  the  case  of 

p  pi  pe 

interest . 

The  dispersion  relationship  (15)  D((i),k)  =  BE  -  C  -  0,  in  the  limit 

k  '  1/r^  >>  k,  is  characterized  Just  by  two  dimensionless  parameters: 

2  2 

8  a  V  /c  and  n  a  (u.,,/u)  .  The  first  parameter  reflects  the  beam  energy  (Ek) 
o  bll  p  D 

and  n  the  density  ratio  of  the  beam  and  the  background  plasma.  This  can  be 
seen  by  normalizing  all  frequencies  with  respect  to  and  all  wave  numbers 
with  respect  to  Wp/v^.  It  is  for  this  reason  that  we  are  able  to  classify  the 
beam  plasma  instability  according  to  just  two  parameters  8  and  n.  over  a  wide 
range  as  discussed  below. 

E- tensive  numerical  and  analytical  studies  of  the  dispersion  relationship 
(15)  lead  to  the  classification  shown  in  Fig.  4.  In  this  figure,  the  n  " 
plane  is  divided  into  five  domains.  In  each  of  these  domains,  the  growth 
rates  of  the  Instability  scale  differently  but  the  instability  mechanism  can 


be  described  simply.  For  example,  the  hybrid  mode  examined  in  the  previous 


section  falls  in  Domain  V  of  Fig.  4.  We  shall  now  document  the  approximate 
growth  rates  for  these  domains.  We  should  also  point  out  that  the  boundaries 
between  the  domains  are  not  sharp.  The  main  use  of  this  classification  is 
then  to  provide  an  immediate  determination  Df  the  Instability  of  a  given 


system,  once  the  beam  energy  and  the  Dacxground  plasma  d'^nsitles  are  given. 


Some  examples  will  be  given  in  the  next  section.  In  all  these  examples,  the 


approximate  formulas  given  below  have  been  favorably  compared  with  the  full 
dispersion  relationship  (15). 


(A)  BY  <  1  (i.e.,  <  212  keV). 


If  BY  <  1,  the  beam  is  essentially  non-relativistic.  In  this  case, 

C  «  0  in  (15),  and  the  dispersion  relationship  for  the  classical  two-stream 
instability 


1 


(19) 


is  readily  recovered  from  (15),  regardless  of  the  value  of  k  .  We  subdivide 

2  2 

into  two  cases:  n  “  oj.  ,,/6j  <1  and  n  >  1. 

bll  p 


(A1)  n  <  1 


(DOMAIN  I) 


In  this  case,  the  background  plasma  density  much  exceeds  the  beam 
density.  This  is  the  typical  case  treated  in  most  text  books.  The  dispersion 
diagram  is  shown  in  Fig.  5.  According  to  (19),  the  peak  growth  rate  u)j^(max), 
the  cutoff  wave  number  k^,  and  are  given  by 


(0.  (max)  -  u  (/3/2)  (n/2)’^^  -  (20a) 

ip  p 

<3/2)  (20b) 


This  instability  is  absent  for  k  >  k^.  The  case  (A1)  is  labeled  as  I  in  Fig. 
4;  the  dominant  instability  being  the  weak  beam  classical  two  stream 
instability. 

(A2)  n  >  1.  (DOMAIN  II) 

In  this  case,  the  background  plasma  has  a  lower  density  than  the  beam. 
The  dispersion  relationship  (19)  may  also  be  analyzed  using  the  standard 
technique.  The  dispersion  diagram  is  sketched  in  Fig.  6.  The  peak  growth 
rate,  tOj^(max),  the  cutoff  wave  number  k^,  and  for  the  unstable  mode  are: 


(max)  -  0.687  oj 

P 

(21a) 

-  (Up/v^)Cp^^^  +  (3/2)p’^^] 

(21b) 

-  (kv^)/n. 

(21c) 

The  scalings  given  in  (21)  are  quite  different  from  those  in  (20).  Note 
however,  that  (20)  and  (21)  yield  the  same  result  if  p  -  1.  Equations  (21) 
characterize  the  parameter  space  labeled  II  in  Fig.  4,  where  the  beam 
dominated  classical  two  stream  instability  is  prevalent. 

The  dispersion  relationship  (19)  is  well-known.  It  is  unnecessary  to 
furnish  here  a  simple  physical  picture  to  support  the  derivation  of  (19).  We 
now  turn  to  the  case  where  the  relativistic  effects  are  important. 

(B)  BY  >  1  (i.e.,  E.  >  212  keV). 

b 

When  the  beam  energy  is  sufficiently  high,  the  relativistic  effects,  the 
electromagnetic  effects,  and  parallel  and  transverse  motions  are  all  coupled. 


2  2 

For  BY  >  1,  we  may  subdivide  into  three  domains  according  to  n  <  1/6  Y  , 
1/B^Y^  <  n  <  6^Y^,  and  n  >  6^Y^.  These  three  domains  are  labeled  III,  V,  IV, 
respectively.  The  approximate  dispersion  relationships,  together  with  the 
growth  rates  and  k^,,  for  each  of  these  domains  are  displayed  as  follows. 

(B1)  n  <  1/6^Y^  (DOMAIN  III) 


If  we  normalize  ui  -  u/Up,  P  -  JJ/Upt  Eq.  (22a)  may  be  written  as 


2  2 

1  6  Y^'n 

1  -  Z2  “  ^2^ 


(22b) 


2  2 

where  the  parameter  6  Y  n  enters.  Recall  that  we  have  in  this  subsection 
2  2 

assumed  B  Y  n  <  1 . 

The  dispersion  curve  according  to  (22)  is  sketched  in  Fig.  7.  The 
maximum  growth  rate  ^^(max),  the  critical  k^,,  and  are  given  by 
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(23a) 


u^(max)  -  0.687 


%  *  ^  (3/2)(sVn)^^^] 


“r.  * 

r  o 


(23b) 

(23c) 


In  (23a)  and  (23b), 


2  2 
k  G 


n[ 


2  2  2  2 

k  c  +  <»> 

1  P  bj_ 


]. 


(24) 


2  2  2  2 

which  reduces  to  n  In  the  limit  k  c  >>  w  +  o).  .  Note  from  Fig.  7  that  as 

j_  p  bj_ 


k  ->■  0,  +  67i/f}(i)p,  which  is  non-zero.  This  mode  has  been  called  a  Weibel 


1  O  1  ii 

mode,  and  is  predominant  in  Domain  III  in  the  classif ication  shown  in 


Fig.  4. 


(B2) 


2  2 
n  >  3  r 


(DOMAIN  IV) 


This  limit  is  the  opposite  as  the  case  (B1).  [cf.  replace  n  by  1/n]. 


this  limit  we  may  approximate  B  »  1  -  u.  E  ■  -  (k^c^  +  u>^  +  m  ^)/a)^. 

Dll  j_  p  b_L 

C  -  -  and  (15)  becomes 


(1) 


1  - 


bll 


2  2 
k  c 

1 


•22  2 

k  c  -^0)  +0) 

1  P  b_L 


2  2  2 
B  u)  u) 

_]  _P_^ 

2  2  2  ’ 
CO  0 


(25) 


2  2  2  2 

For  simplicity,  we  further  assume  k  c  >>  u  +  ux  to  reduce  (25)  to  an  even 

1  P 


simpler  form  which  is  independent  of  k  ; 


The  solution  to  (26)  is  sketched  in  Fig.  8;  the  maximum  growth  rate 


ui^(raax),  k^,  and  are  given  by 

0),  (max)  -  0.687  oi  (6T)^^^n^^^  (27a) 

i  P 

k^  -  (Vo/ii-p)Cn’^^  +  (3/2)(3Y)^^^n^^^]  (27b) 

u  «  4  (27g) 

r  2  o 


The  dispersion  relationship  (26)  appears  to  be  new.  To  gain  some 
understanding,  we  derive  it  simply  as  follows.  Upon  dropping  in  Eq.  (4), 
we  obtain 


Ym  5 
e^e 


Eli] 


or  equivalently 

-  ‘  ^ll> 

where  we  have  used  the  fact  that  the  self  field  E^g  is  canceled  by  the  Lorentz 

2 

force  (Vq/c)  Sle  by  the  factor  1/Y  .  [cf.  Fig.  1].  Again,  neglecting  E^j^  in 
comparison  to  E^^g  in  (5),  we  have 

-  (28b) 

We  now  use  the  expressions  E,  -  4iT|e(n.E  and  E, ,  =  -  4Triein  E.  in  (28a)  and 

le  ''be  li  ''pi 

(28b).  The  dispersion  relation 
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1 


2 


(29) 


2  2 


0)  0) 


2 

is  then  readily  obtained  from  (28a),  (28b).  Aside  from  the  factor  g  (which 
has  been  assumed  to  be  closed  to  unity  in  this  subsection),  Eq.  (29)  is 
identical  with  (26).  This  simple  derivation  may  be  used  to  provide  a  physical 
explanation  of  the  instability,  similar  to  that  given  in  Fig.  3.  Note  that 
while  this  derivation  is  very  similar  to  that  given  for  the  hybrid  mode,  the 
scalings  of  the  growth  rates  for  the  two  instabilities  are  different  [cf. 

Figs.  2,  8;  Eqs.  (13).  (26)]. 

(B3)  <  n  <  (DOMAIN  V) 

2  2  2  2 

The  remaining  domain  is  l/g  Y  <  n  <  6  Y  .  This  is  labeled  as  Domain  V 

in  Fig.  4.  The  dominant  instability  for  this  parameter  space  is  the  hybrid 

mode  studied  in  detail  in  Section  II.  The  simple  dispersion  relation  was 

given  by  Eq.  (13)  and  sketched  in  Fig.  2  of  that  section.  We  note,  from  the 

simple  derivation  given  in  Sec.  II,  that  this  domain  is  characterized  by 

Y  >>  1  and  >>  1  [cf.  Eqs.  (1)  and  (14)].  The  condition  >>  1  is 

Dl  P  Dl  P 

just  n  >>  1/6^Y^  when  g  -»•  1 . 

Finally,  we  note  that  the  classification  shown  in  Fig.  4  is  symmetrical 
about  n  -  1 .  This  is  reasonable  if  one  views  the  beam  plasma  interaction  in  a 
frame  moving  at  the  mean  speed  of  the  two  species.  We  also  stress  that  the 
simple  dispersion  relations  given  for  Domains  (I)  -  (V)  are  strongly  supported 
by  the  numerical  computation  of  the  full  dispersion  relationship  (15). 

IV.  Some  Examples 

In  this  section,  we  give  a  few  examples  on  the  use  of  the  classification 
schemes.  These  examples  are  drawn  from  laser  accelerations  of  ultra-high 
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energy  particles, the  plasma  assisted  injection  and  extraction  of  electrons 
into  a  modified  betatron.^  The  simple  formula  and  Fig.  4  provide  ready 
characterization  of  the  instability,  in  terms  of  the  growth  rate  and  the 
critical  axial  wave  number  (above  which  the  system  is  stable).  One  example  is 
given  to  each  of  the  Domains  III,  IV,  V.  [Domains  I  and  II  cover  the  more 
familiar  plasma  interactions  with  a  non-relativistic  electron  beam,  and  will 
not  be  considered  further  in  this  section]. 

2 

A.  Laser  Beat  Wave  Acceleration 


Consider  the  acceleration  of  an  electron  beam  whose  energy  is  10  MeV. 

Let  rjj  =  beam  radius  -  0.5  mm.,  Ij^  =  beam  current  -  0.1  A.  This  beam  is 
embedded  in  a  background  hydrogen  plasma,  of  density  np  =  10^ ^/cm^, 

7  1 2 

say.-^  These  parameters  yield  Y  -  20.6,  u  =  w  -  5.64  *  10  rad/sec,, 

7  ft  P  ? 

1  ■  3*11  *  10  rad/sec,  u.  •>  6.4  *  10  rad/sec.  Then  n  »  oj.  , , /oj 
bii  bj_  bll  p 

-  3.03  X  10  which  is  much  less  than  1/Y^B^  =»  2.36  *  10  Thus  the 
parameters  lie  within  domain  III. 

The  growth  rates  and  critical  wave  numbers  may  be  evaluated  using  Eqs. 

(23).  We  arbitrarily  take  k  •  1/r.  .  Then  (24)  yields  n  “  n  k  c/m 

±6  IP 

»  0.106n  ■  3.22  X  10"^^.  Thus,  Eq.  (23)  yields  -  2.03  >«  10^  sec"’.  For  a 

pulse  length  Tp  -  3-33  »  10  sec,  “jfp  -  0.68,  which  is  the  total  number  of 

e-folds  during  the  beam  pulse.  Thus,  the  instability  is  very  mild.  The  half 

wavelength  at  k„  is  A  /2  -  tr/k  -  0.017  cm,  according  to  (23b). 

c  c 

The  instability  growth  rates  [cf.  (23a)]  is  proportional  to  the  1/3  power 
of  Ijj/Y.  Thus  the  conclusion  given  in  the  previous  paragraph  regarding  the 
unimportance  of  the  instability  remains  valid  for  a  substantially  higher  beam 


current,  or  at  a  higher  beam  energy 


B.  Plasma  Assisted  Injection  of  Electron  Beam  in  a  Modified  Betatron 


Configuration 

The  modified  betatron  configuration^  consists  not  only  of  the  usual 
betatron  magnetic  field  but  also  a  strong  toroidal  magnetic  field.  To  inject 
an  electron  beam  into  the  modified  betatron,  it  is  proposed  to  create  a  plasma 
channel^  into  which  the  electron  beam  propagates.  In  the  ion  focus  regime, 


0)  •  0)  . 

P  Pi. 


To  determine  the  extent  of  the  beam  ion  interaction,  consider,  for 


example,  an  electron  beam  of  energy  1  MeV,  of  beam  current  1  fCA  and  beam 


.11 


radius  of  0.5  cm.  Then  6  =  0.9^1,  Y  =  2.957,  n^^  •  2.7  x  10  /c.c.. 


9  1 0 

=  5.7  X  10  rad/sec.,  =  1.69  x  10  rad/sec.  If  the  ion  channel 


“bli 


is  formed  by  benzene  gas,  m,  -  73  x  18^10  m  .  Take  n.  »  n  ,  then 

“  ’  1  e  1  e’ 


7  2  2'^ 

u  =  0)  ,  *  7.4  X  10  rad/sec.  Thus  n  =  oj.  ,,/a)  »  5-8  x  10  .  Together  with 

p  pi  bll  p 


=■  1  MeV,  we  note  that  the  parameters  lie  within  Domain  IV  [cf.  Fig.  4],  If 


2  2  2  2 

we  arbitrarily  take  k  -  1/r.  ,  then  k  c  >  oi  +  to.  and  the  dispersion 

1  1  P 


relation  (26)  is  a  good  approximation  to  (25). 


The  growth  rate  and  k^,  are  given  by  Eqs.  (27a)  and  (27b).  Using  the 


8  —1 

parameters  given  in  the  previous  paragraph,  we  deduce  =  4.3  x  io  sec”  ,  so 


that  during  the  time  scale  tp  *  20  ns,  say,  corresponding  to  a  single  turn 


1 

injection  time  of  the  NRL  modified  betatron,  the  total  number  of  e-folds  of 


the  instability  is  ”8.6.  This  suggests  that  the  beam  type  Weibel  mode 


may  be  marginally  important  since  the  above  analysis  excludes  stabilizing 
influence  such  as  betatron  oscillation,  external  magnetic  field,  random  ion 
motion,  '  the  enharmonic  property  of  the  channel,  and  the  convective  nature  of 
the  mode.  Thus,  the  total  number  of  e-folds  is  expected  to  be  somewhat  less 


than  8.6,  and,  as  a  result,  be  tolerable.  From  (27b),  we  have  X^/2  =  ir/kj. 


13*4  cm.  This  provides  an  estimate  of  the  axial  extent  of  the  mode. 


C.  Plasma  Assisted  Extraction  of  Electron  Beam  In  a  Modified  Betatron 
Configuration 

The  idea  of  using  a  plasma  channel  may  also  be  adopted  for  the  extraction 
of  an  electron  beam,  after  the  latter  has  completed  the  acceleration  phase  in 
a  modified  betatron  configuration.^  Note  that  such  an  idea  is  especially 
attractive  for  beam  extraction,  because  the  contamination  of  the  chamber  would 
not  become  an  issue . 

To  determine  the  interaction  between  the  post-accelerated  beam  and  the 
channel  ions,  consider,  again,  an  electron  beam  of  current  1  KA,  beam  radius 
Tjj  =  0.5  cm,  this  again  yields  n^^  =  2.7  X  lO^Vc.c.  If  benzene  is  used,  and 

7 

if  n.  »>  n^,  then  uj  =a).“7.4*10  rad/sec  as  in  Subsection  IV  B.  If  the 
lb’  p  pi 

7 

energy  of  the  beam  is  about  50  MeV,  say,  (Y  *  100),  then  ^ 

9  2  2 

rad/sec  and  oj.  -  2.9  *  10  rad/sec.  The  parameter  n  =■  -  0.152. 

bjL  bll  p 

These  values  of  (E^,  n)  lie  within  Domain  V  in  Fig.  4.  Thus,  the  post- 

accelerated  beam  is  subject  to  the  hybrid  mode. 

The  growth  rate  and  for  the  hybrid  mode  are  given  by  (11).  Using  the 

8  “1 

above  parameters,  we  have  o),  =»  oj  lo,.  “  4.6  *  10°  sec  .  This  growth  rate  is 

i  p  bj_ 

similar  to  that  in  Subsection  B.  However,  here  X  /2  •  n/k  »  1m,  which  is 

c  c 

quite  long.  In  fact,  for  parameters  similar  to  the  NRL  modified  betatron,  the 
length  of  the  channel  is  always  less  than  60  cm,  by  geometrical 
considerations.  The  relatively  mild  growth  and,  more  importantly,  the  long 
wavelength  nature  of  the  hybrid  mode,  render  this  method  of  beam  extraction 
attractive. 

V.  Discussiona 

In  this  paper,  we  classify  beam  plasma  interaction  according  to  just  two 
parameters:  the  beam  energy  and  the  beam-plasma  density  ratio.  The 
dispersion  relationship  Is  fully  electromagnetic  and  fully  relativistic,  and 
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has  been  given  by  various  authors.  Here,  we  divide  the  parameter  space  into 
several  domains,  in  each  of  which  the  instability  growth  rate  and  axial  scale 
lengths  are  simply  given.  We  stress  the  simplicity,  the  exposition  of  the 
instability  mechanism,  and  the  ease  with  which  such  a  classification  may  be 
applied  to  a  particluar  situation.  Examples  are  given. 

While  the  classification  given  in  this  paper  provides  an  immediate 

assessment  of  the  instability  growth  rates  -  over  a  wide  range  of  beam 

energy  and  of  number  densities,  the  following  limitations  should  be  kept  in 

mind,  (a)  The  role  of  external  magnetic  fields  have  been  ignored.  In  the  ion 

focus  regime,  this  assumption  may  be  argued  to  be  valid.  However,  in  many 

other  cases,  this  assumption  cannot  be  made  and  the  classification  given  here 

can  no  longer  be  used.  Intuitively,  the  magnetic  field  may  be  ignored  if 

2  2 

n.  <<  01.,  where  Q  and  oi.  are  respectively  the  relativistic  cyclotron 
d  1  J  d 

frequency  and  the  relativistic  plasma  frequency  of  the  species,  (b) 
Perhaps  equally  restrictive  is  our  over-simplification  of  the  equilibrium 
fields.  The  influence  of  the  self  fields  on  this  classification  remains  to  be 
determined,  (c)  When  the  densities  are  high,  collisional  effects,  which  have 
been  neglected  in  this  paper,  may  become  important,  (d)  We  have  ignored 
betatron  oscillations  of  the  beam  electrons,  and  possible  ionic  motions  in  the 
unperturbed  states.  These  effects  are  likely  to  be  stabilizing,  (e)  We  have 
not  characterized  the  absolute  or  convective  behavior  of  the  instabilities. 
From  previous  experience  with  beam  instabilities,  we  anticipate  that  the 
Instability  is  mildly  absolute,  whose  time  asymptotic  response  behaves  like 
exp  (t^),  where  o  <  a  <  1,  typically.  However,  inclusion  of  betatron 
oscillations  would  tend  to  make  the  instabilities  convective.^  (f)  We 
conjecture  that  the  details  of  the  beam  profiles  are  not  very  crucial  in  the 
classification  of  beam  instabilities.  Experience  suggests  that  the  radial 
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scale  length  of  the  beam  provides  a  convenient  estimate  for  k~^ .  The  details 

1 

would  be  less  imporcant.^”^ ^  Incidentally,  in  the  classical  two  stream 
instability,  the  presence  (or  absence)  of  discrete  modes  depends  on  the 
sharpness  of  the  boundaries.  However,  the  time  asymptotic  response  is  just 
equal  to  that  associated  with  the  highest  local  growth,  regardless  of  the 
profile.^*^^ 

Finally,  we  comment  on  a  few  works  on  electron-ion  instabilities. 
Manheimer^^  provided  a  detailed  calculation  with  the  use  of  a  self-consistent 
equilibrium  distribution  for  the  ions.  However,  the  electron  inertia,  which 
is  extremely  important  for  the  hybrid  mode,  is  neglected  (justifiably  because 
of  the  strong  magnetic  field  he  assumed).  Buchanan^®  includes  the  electron 
inertia,  but  apparently  he  ignored  the  possible  cancellation,  between  the 
and  (Vg/c)Bj^g,  to  order  1/Y  ,  as  explained  in  Sec.  II  and  Subsection  Illb. 
Davidson  and  Uhm^^  included  the  case  of  a  strong  magnetic  field,  and  Uhm^° 
recently  concluded  that  two  stream  instabilities  are  unimportant  for  highly 
relativistic  beams.  Godfrey^^  and  his  coworkers  performed  particle 
simulations  in  the  ion  focus  regime.  None  of  these  authors  gave  the  hybrid 
mode  nor  the  classification,^^  as  stressed  in  the  present  paper. 
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Eb  (keV) 


5.  4  —  Classification  of  beam  plasma  interaction.  On  the  boundary  between 
Domain  IV  and  V,  tj  *  and  on  that  between  III  and  V,  1/tj  —  b\^. 
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le  normalized  (complex)  frequency  as  a  function  of  the  normalized 
axial  wave  number,  for  Domain  1  of  Fig.  4. 
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Appendix 

HYBRID  MODE  IN  A  FINITE  BEAM 


In  this  appendix,  we  consider  in  some  detail  the  Interaction  between  a 
relativistic  electron  beam,  of  radius  r^,  and  a  background  plasma,  the  latter 
being  infinite  in  extent.  We  focus  on  the  hybrid  mode.  We  show  that  the 
growth  rate  given  by  Eq.  (13)  i3  also  applicable  when  the  beam  is  finite. 

This  study  provides  additional  confidence  in  the  replacement  of 

by  p/r^j,  where  p  is  some  quantized  number  of  order  unity,  or  higher. 

When  the  equilibrium  self  fields  are  ignored,  the  interactions  of  the 
electron  beam  and  a  cold  background  plasma  are  described  by  the  following 
coupled  differential  equations: 


i  {Hi  -  ifig  -eE  -  — 

r  dr  ‘  6  dr  '  2  6  z  ^  z  2 

r  r< 


6  dr 


z  ^ 
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Here,  E^  and  B^  are  the  rf  field  components  of  the  TM  and  TE  modes,  which  are 
coupled  if  the  azimuthal  mode  number  H  •  0.  Small  signal  dependence  of  the 
form  f(r)  exp  (iut  -  lie  -  ikz)  has  been  assumed.  The  symbols  in  (A1 )  and 
(A2)  are  defined  as 
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2  2 

kcoj  (jo  (kc  -  u)V  /c) 

.  -L  [_p  .  _ !_]. 

6c 


(A5) 


Consider  first  a  local  analysis  of  (A1).  Take  £  -  0  for  convenience  and 
replace  d/dr  by  -  ik^.  Then  Eq.  (A1)  yields 


k  ;/6  -  e  -  0 

1 


(A6) 


which  can  easily  be  shown  to  be  the  dispersion  relationship  (15).  It  is  also 
quite  clear  that,  in  general,  l/rj^  provides  a  natural  transverse  scale  length 
of  the  perturbations.  In  fact,  one  may  use  a  normalized  variable  x  =  r/r^  to 
render  (A1)  and  (A2)  dimensionless. 

For  general  £,  and  for  a  constant  beam  density  up  to  r  -  (and  zero 
for  r  >  rjj),  a  transcendental  dispersion  relationship,^^  in  terms  of  Bessel 
functions,  may  be  obtained  from  (A1 )  and  (A2).  This  dispersion  relationship 
may  be  obtained  by  requiring  that  and  be  continuous  at  r  »  rj^,  but 

f  f 

E  and  B  be  discontinuous  by  an  amount  which  is  to  be  computed  from  an 
z  z 

integration  of  (A1 )  and  (A2),  due  to  a  discontinuity  in  p.  This  dispersion 
relationship  may  be  simplified  considerably  under  the  assumption 
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<< 


In  this  case,  it  reads 
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where 
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(A9) 
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and  X  is  the  v-th  zero  of  the  Bessel  function  J. .  To  the  lowest  order  in 

XV  i 

iDj^^Pj^/c  [cf.  (A7)],  (A9)  then  yields 


2  2  2 
u)  S  u) 

1  -  -i— U.  -  0. 
2  2 
n  01 


(A10) 


The  solution  to  (A10)  is 


i  (All) 

Which  is  just  the  dispersion  relation  (11)  for  the  hybrid  mode  [aside  from  the 
factor  6,  which  has  been  assumed  to  be  closed  to  unity  in  the  derivation  of 
(11)].  The  dispersion  relationship  (All)  is  sketched  in  Fig.  2. 

Finally,  we  remark  that  as  long  as  (A7)  is  valid,  the  growth  rate  given 
by  the  local  analysis  is  anticipated  to  be  valid  also  for  a  diffuse  beam 
profile.  This  has  also  been  pointed  out  by  Tajima.^® 
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